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> , 

Q , Decay law of a complicated unstable state formed in a high energy collision is described by the 

!Z ■ 

^^ ' Fourier transform K{t) of the two-point correlation function of the scattering matrix. Although 

each constituent resonance state decays exponentially the decay of a state composed of a large 
number A^ ^ 1 of such interfering resonances is not, generally, exponential. We introduce the decay 
r~| I rates distribution function w(r) by representing the decay law in the form of the mean- weighted 

^ ' decay exponent K{t) = f^ dV e w{r). In the framework of the random matrix approach we 

investigate the properties of the distribution function w{T) and its relation to the more conventional 
statistics of the decay widths. The latter is not in fact conclusive as concerns the evolution at the 
times shorter than the characteristic Heisenberg time. Exact analytical consideration is presented 
for the case of systems without time reversal symmetry. 

PACS numbers: 05.45.Mt, 03.65.Nk, 05.45. +b, 24.60.-v 
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I. INTRODUCTION 

The temporal aspect of chaotic resonance scattering is an interesting and important issue 
which repeatedly attracted much attention starting from the seminal works of Wigner and 
Smith l|, |2|. Later on different problems concerning the duration of resonance collisions 
were posed and discussed in detail. Distribution of resonance widths on the one hand and 
statistics of partial and proper dela y ti mes on the other have been investigated in a number 



statistics or partial andproper dela y ti mes on tne otne 
of papers [a, U, U U U U U M lul, Id, U MUM. 



It is generally known that a single quasi-stationary state formed by the time t = in 
a collision with a given energy E decays afterwards exponentially with the decay rate F^ 
which is given by the imaginary part of the resonance complex energy Sr = E^ — ^Tr- In 
accordance with the Bohr energy-time uncertainty principle this rate defines the width of 
the Breit- Wigner resonant curve. 

However at high enough collision energy E a very complicated unstable state is formed 
which is a superposition of many, A^ ^ 1, interfering resonance contributions whose spectrum 
is given by the N complex eigenvalues Sr = Ej. — ^F^ oi N x N matrix of the non-Hermitian 
effective Hamiltonian Ti = H — ^W, W = AA^ . The resonances, generally, overlap so 
that their mean width is larger than the mean level spacing, (F) > A. The evolution of such 
an unstable configuration formed via an incoming channel a and decaying onto an outgoing 
channel b is described [4, l7| by the Fourier transform 

1 f°° 
K^«(r) = — / duje'^^C'"'{uj) (1) 

of the two-point S-matrix correlation function 

Here T'^ = 1 — \S'^'^\'^ is the transmission coefficient in the channel c and tn = 27r/A is 
the Heisenberg (Weisskopf) time. It is convenient to measure the time in the units of 
this characteristic time interval, r = t/tn- Correspondingly, the dimensionless energy shift 
uj = etH is introduced in the Eqs. ([H [2]). To avoid complications irrelevant in the present 
context we restrict for a while our consideration to inelastic collisions only. 

Contrary to the case of an isolated quasi-stationary state, the function K^°'{t) does not, 
generally, decays exponentially and a connection to the widths of the constituting resonances 



is not immediately seen. More than that, an opinion has not once been expressed that the 
very notion of the resonance widths becomes irrelevant when resonances strongly overlap. 
Nevertheless, intuitively, one would expect that decay properties of a complicated unstable 
state should somehow depend on the distribution of the individual resonance poles of the 
scattering amplitudes in the complex energy plane. Could then a non-exponential decay 
results from an incoherent mixture of exponential decays with different rates? If and how 
such a mixture is related to the statistics of the widths of resonances? We address these 
questions below. They are answered explicitly by the example of the chaotic systems without 
time-reversal symmetry. A rigorous analytical solution is found in this case. 

II. DECAY RATES DISTRIBUTION FUNCTION 

The function C'^°'{uj) is analytical in the lower part of the complex u plane and, obviously, 
satisfies the condition C^°-{—uj) = C^"" {uj) . As a result its Fourier transform K^^-ij) = 
when r < and is real when r > 0. More than that, independently of the number M of the 



reaction channels, 17 1 



K'"'{0) = 1 and K'"'{t) > when r > . (3) 

Notice at last that the mean reaction cross section is expressed as 

((T^") = T^ T" C^"(0) =T^T'' / dr K^\t) . (4) 

Jo 

Such properties suggest the Laplace representation for the decay law function 

/■oo 

K'"'{t) = / ci7M;''"(7)e-^". (5) 



^0 

(In the ordinary units 7 = Ft/^ = ^ F where F is the decay rate.) In view of the Eq. ([3]) 
the function w'^"' is normalized to unity. 



rf7w''«(7) = J^f"^(0) = l. (6) 







It follows from the Eqs. ([T], [5]) that 



c'«(.).r,,^ (7) 

Jo ^ + IU} 



so that the correlation function is not, generally, single-valued in the complex uo plane. 
Finally, according to the Eqs. (jH [7]), 

(a''") = T'' T" r ^7^^^^ . (8) 

Jo 7 

From this point on we suppose that all channels are statistically equivalent, T^ = T"" = T, 
so that the functions C(u;), K{t), w{'~f) do not depend on the channel indices. 



III. THE SEMICLASSICAL ASYMPTOTIC EXPANSION 

Let us consider first the semiclassical limit: A^ ^ 1, M ^ 1 but the ratio m = M/N , is 

finite though small, ?n <C 1, (this condition is physically justified). The function ^'(7) can 

easily be extracted from the results reported in: J5|, M\ (GOE) 

1 7(^) 1 2T MT^ 

MT 7(c<;) + iuj 'yw + i^ {"^w + i^Y {iw + i^Y 

where 7(0^) is a slow varying function and 7(0) = 714/ = MT = tu^w where Tw = ^ MT is 

the Weisskopf width. The asymptotic series in the r.h.s which is obtained by expanding near 



the pole uj = i'jw coincides with that derived first in different manner in: [18|. The first term 

corresponds to the Hauser-Feshbach approximation and implies the purely exponential decay 

Q-iwr Tf^iiQ-j^ tiie subsequent ones account for deviations. In the limit TV, M^oo, 0<rra<^l 

the Weisskopf width equals to the empty gap between the real energy axis and the cloud 

of the poles corresponding to contributing resonances |9|. The width Tw is in this case the 

smallest possible width. 

It immediately follows from the expansion (Q that 

MT'^ 
w{^) = 5(7 - 7iy) - 2T<5'(7 - 7w) + ^- S"{i - ^w) + - • (10) 

The function ^(7) thus obtained is neither smooth nor positive definite. We will see below 
that the expansion iQ is in this respect somewhat misleading. Nevertheless it works well in 
certain cases. In particular, for the mean cross section we get from the Eqs. ([HI flOjl 



<->'°°^' - F (1 - i^ ^ ■■■) ■ (") 

Similar calculation yields in the GUE-case {cr)^'^^^^ = jj . 

Notice that the difference (o")^*-^^^^ — {a)^'^'~'^^ = -^ determines the "weak localization" 
in the quantum transport. 



IV. DECAY RATES VERSUS WIDTHS STATISTICS 



Heuristic" arguments have been adduced in [19| which assume a simple connection be- 
tween the function w{'~f) and the distribution Pm{i) of the resonance widths 

WMh) = A'y^ PMh) (12) 

with some normahzation constant A. A formula of such a kind has been derived a little later 
in 8|. It was shown, however, that the relation supposed can reproduce only the long-time 
asymptotic behavior. Only the most long-lived resonances survive by this time, so they can 
be expected to obey the Xm widths statistics. 

Pm(7) = -OT)!^ (/5^)'^e-^-, (^) = MT = ^w, P = l,2. (13) 

Here the parameter (3 marks the Dyson's symmetry class: GOE (/3 = 1) or GUE {[3 = 2). 
Notice that in our dimensionless units r/(r) = '-y/MT which means that (7) = MT = jw 
In the case of isolated resonances the Weisskopf width defines the mean rather than the 
minimal width. The formal extending the relation flT2|) to all values of 7 yields the result 



KM=Am(m + ^] ^\, ^^^ (14) 

which reproduces rightly the power law asymptotic behavior |5| though the proper calcula- 
tion of the constant is beyond the validity of this approximation. The constant A can be 
chosen to quantitatively fit the correct asymptotics. However at shorter times the found in 
such a way expression disagree with the actual decay law, the discrepancy being the stronger 
the closer is the transmission coefficient T to its maximal value 1 (see |8|). There is no way 
to reconcile all necessary conditions which must me satisfied by the decay function Km{t) by 
means of only one matching parameter A. There are two possible reasons why the supposed 
relation flT^ fails: i. the x\i distribution is not valid when resonances overlap; ii. the naive 
derivation is wrong. We will clarify below these guess-work by the example of the systems 
with broken time-reversal symmetry which are described by the Gaussian unitary ensemble 
(GUE) of random Hamiltonians. 



V. SYSTEMS WITH NO TIME-REVERSAL SYMMETRY 



General consideration 



The GUE two-point correlation function derived in: 20| reduces in the case of equivalent 
channels to 



^y^) - Jo "-^Jo "^1^ A+Ai (l+TAi)«+2 ^ 

{ (I - 1) (1 - rA)(l + TAi) - i-^ [(1 - T\) + (1 + TAi)]} 



(15) 



Making use of the following simple identity 



oo 



ATTA-y, d,e-*.-) (16) 

we factorize the integrations over A and Ai. On the next step we use a second identity 

The correlation function (fT5!) is thus reduced to a linear combination of terms like 



-j /"CXD /"OO /■! 



-i(jjX 



After a chain of simple transformations of the integration variables the Fourier transform of 
the latter function reduces finally to 

The symbol G stands for the step function. All successive integrations in this expression 



can be carried out explicitly 2ll ]. 



To simplify further consideration we restrict our calculation to the case of the perfect 
coupling to the continuum, T = 1, when the naive consideration is the least satisfactory. 
Only the term with fi = fi' = remains in this limit and 

-| i>oo POO /■! ~iu)X 

CMico) = J7-, dj Hi"' e-^' / d\{l- \r-' e-^' ^t"- • (20) 

Ml Jo Jo Jo 7 + 7 + «^ 

The subscript M explicitly indicates the number of open channels. Notice that when T = 1 

the formula (l20l) and so all consequent results are valid also for the elastic collisions including 

the purely elastic process with one open channel M = 1. 
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The found representation enables us to evaluate the mean cross section in the following 
simple and elegant way: 

1 /"OO 1 POO 1 



{a) = CMiO) = — d^FM{l)^M^iil) = j^ rf7e-^Fo(7) = — . (21) 

The third equality follows from the fact that the functions 

POO -^T) p 

Fm{i)= dvv''^. ^M-i{l)= da^-CV'~'e-< (22) 

Jo 1 + ^ Jo 

which appear in this equation obey the simple recursions 

Fm{i) = -^FM^ih), :^$A/(7) = M$M_i(7)- (23) 

(27 07 

The decay function looks now as 

/oo in n' 

c?7e-^^^y di j duQ{r-l + u/^')u''-'e-\ (24) 

The later consideration depends on whether the time r exceeds the Heisenberg time r^ = 1 
or not. 

B. Long-time asymptotics 

After the Heisnberg time, r > 1 , the step function in the Eq. (l24l) equals to one in the 
whole integration region and the decay function gets the required form 



oo 



(26) 



Km{t>1)= rf7e~^"w;>(M,7) (25) 

Jo 

with the decay rates distribution function given by 

«7>(M,7) = ^/;rf7VM-i(7'), 
ipM-ib) = fjduiy'^-^e-'' = (M- 1)! -r(M,7) 

where r(M, 7) is the incomplete F function 

r(M,7) = (M-l)!e-^^ 1-. (27) 

■^^ m! 

m=0 

The next integration over 7' yields finally 

w,iM, 7) = r(M+i..)-.r(M,.) + ^^ _ 1 = ^^^(^) _ (1 _ ^^) , 
^?^/(7)=e--Eito (l " f ) S = e--P(M-i)(7) • 



(28) 



As expected, this distribution is smooth and positive definite. Notice that, due to subtraction 
of the last two terms, the distribution tf;>(M, 7) vanishes as 7^ when 7 —> 0. 
At last the final integration over 7 yields 



thus fixing not only the power but also the constant of asymptotic power behavior. 

The found results are closely connected to statistics of the decay width. Indeed the widths 
distribution function in the case of M equivalent channels, arbitrary degree of resonance 



overlapping and perfect coupling to the continuum is currently well known to be 11 1 

1 n 



and has nothing to do with the XM-distribution (IT3l) It is immediately seen that 

^M-iil) = (M - 2)\^^pM-ih) (31) 

(see the Eq. f l26|) ). This results in the following connection between the two distributions 



(32) 



1^ b'pMil) + lo' dl' 1" PMil')] =^{i + l) Io'dl'l"pM{l') ■ 
The extra term which appears in this relation is missing in the naive formula (fT2|) . 

C. Short-time behavior 



At the times shorter than th-, r < 1 , an additional term arises, 

rfz/ e(r - 1 + Z//7) u^^-^ e-"" = ^m-i{i) - J dv u^^"^ e'" , (33) 

Jo 

which depends on the time. As a consequence the decay law contains along with the smoothly 

distributed exponential contributions an additional polynomial term 



Km{t < 1) = /o" ^7 e-^^ wm{i) + Ti 



1 l-(l-r)^ 1 



r2 






/J" ^7 e ^" WMin) - if Em=o (-1 . 

m + 2 



Formally, we still can present this formula in the form 

/■oo 

Km{t<1)= rf7e-^^w<(M,7) (35) 

Jo 

with a "distribution density" w<^{M, 7) = WMil) + ^m ^ (7) • The price paid is the singular 

nature of the additional weight function 

-r-h)^-^_i:(-ir[,^^^J^%). M>2. (36) 

This singular "distribution" differs in two important respects from that extracted from the 
Verbaarschot's asymptotic expansion: first, it contains only a finite sum of derivatives of the 
(5-function and, second, all singular terms are concentrated near 7 = rather than 7 = 714/. 
We should stress that the both regular and singular parts play, generally, equally impor- 
tant roles. In particular, 

A^„(0)^^-^^1. (37) 

The first contribution comes from the regular and the second from the singular parts of the 
"distribution" u;<(M, 7). They are of the same order of magnitude when M ^ 1. 
Integration over 7 gives now 

l_ 1 - (1 - T^ ^^' 

M r2(l + r) 
Combining all found results we arrive finally at 



KMir<l) = - '^' , K,M = l. (38) 



^M(-) = ^^,(Y^[l-e(l-r)(l-rY^], 0^r<oo. (39) 

This function satisfies all necessary conditions and is continuous though not analytical in 
the point r = rvi/ = 1- 

In the semiclassical limit M ^ 1 the characteristic decay time tw = l/'^w = 1/M <^ 1 
is much shorter than the Heisenberg time. The polynomial term is therefore of principal 
importance. It becames obvious from the following equivalent presentation of the Eq. ( l38ll 

K^ir < 1) = jv^ [e-A/ln(l+r) _ gA/lnd-r)] 
_Af(,+0(r3)) 2.inh[f(r^+0(r^))] _ 

Thus such an exponential semiclassical decay with the characteristic Weisskopf 's decay rate 
7iy cannot be directly traced to the statistics of resonance widths. Deviation from the 
exponential law due to the neglected terms becomes significant after the time r^ ~ l/yM ^ 



-Af(r+0(r3)) 2.inh[f(;^+0(r^))] _ ,,^ ^ '^^^ (^O) 



Tw |22|,l23|. 



D. Mean cross section 

The mean cross section can be expressed in two equivalent forms 

(^) = 17 = lo dr Knir) + T dr K^ir) = 

r ? ^-' [PiM-D (7) - 1 + ^] + r ^ir -^ b)'-^- 

The long-time contribution (r > 1) rapidly decreases when the number of channels M grows. 
Even for M = 1, 2 it amounts only to 30% and 10% respectively. When M ^ 1 the following 
asymptotic expansion 

1 9 /•Af/2 9 

^ drKMir) = Km{1){^1 dee-(^+^^V2M+...) ^ -i^e^*^'- (42) 

shows that this contribution diminishes very fast. The mean cross section is defined by the 
short time evolution. 

On the other hand, the contributions of the regular and singular parts of the decay rates 
distribution are equally important. Indeed 



, . M 

— e ^ 



P(M-i)(7)-l + ^l=^ + E-- (43) 

The logarithmically growing extra term is perfectly compensated by the contribution of the 
singular part. 

VI. CONCLUSIONS 

The decay law K(t) of a complicated unstable state of a classically chaotic system can be 
formally presented in the form of a Laplace integral as a weighted mean value of the decay 
exponents. By this the notion of the decay rates distribution is introduced irrelative of statis- 
tics of the resonance poles of the scattering amplitudes. The connection of this distribution 
with statistics of the resonance widths is then investigated. Exact analytically solution is 
found in the case of systems with broken time-reversal symmetry. It is demonstrated that 
only the long-time r ^ th = 1 asymptotic behavior of the decay law is governed by the 
statistics of the resonance widths. For the shorter times t < r^ = 1 the contributions of the 
connected to the widths statistics smooth part and the singular part of the " distribution" 
U7<(M, 7) are equally important. In particular, the approximately exponential semiclassical 
decay with the characteristic Weisskopf 's decay rate 7^ results from interrelation of the 
both of them. 
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